#ifndef __RICHALGO_H__
#define __RICHALGO_H__

int richardsalgo(double H0[], int solvecase, int iz, double dz,double dt,
                  double material[][3],double k0[],double phi[],double rhomatrix[][3],int iter_t);

#endif
